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ABSTRACT 



Based on A'-band integral-field spectroscopy, we analyze four Wolf-Rayet stars of the nitrogen sequence (WN) found in the inner part 
of the Quintuplet cluster All WN stars (WR 102d, WR 102i, WR 102hb, and WR 102ea) are of spectral subtype WN9h. One further 
star, LHO 110, is included in the analysis which has been classified as Of/WN? previously but turns out to be most likely a WN9h 
star as well. The Potsdam Wolf-Rayet (PoWR) models for expanding atmospheres are used to derive the fundamental stellar and wind 
parameters. 

The stars turn out to be very luminous, log (L/L©) > 6.0, with relatively low stellar temperatures, T, « 25-35 kK. Their stellar winds 
contain a significant fraction of hydrogen, up to Xu ~ 0.45 (by mass). We discuss the position of the Galactic center WN stars in the 
Hertzsprung-Russell diagram and find that they form a distinct group. In this respect, the Quintuplet WN stars are similar to late-type 
WN stars found in the Arches cluster and elsewhere in the Galaxy. 

Comparison with stellar evolutionary models reveals that the Quintuplet WN stars should have been initially more massive than 
60 Mq. They are about 2.4-3.6 Million years old, and might still be central hydrogen burning objects. The analysis of the spectral 
energy distributions of the program stars results in a mean extinction of A^- = 3. 1 ± 0.5 mag (Ay = 27 ±4 mag) towards the Quintuplet 
cluster. 
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1. Introduction 

The Galactic center (GC) is a unique region to study star forma- 
tion and evolution in the special environment close to the super- 
massive black hole. Three very massive clusters have been found 
within 35 pc projected distance to the GC, hosting a large frac- 
tion of the known Galactic massive stars such as luminous blue 
variables (LBV) and Wolf-Rayet (WR) stars. The Central cluster 
is lo cated within 1 pc from the central black hole ( Krab be et alj 
11991 ). Two further massive star clusters, the Arches and the 
Quintuplet, are located in the vicinity of Sgr A*. Interestingly 
the three clusters diff er in age, the 2.5 Ma old Arches cluster 
(iFiger et a l. 2002, Nai arro et al.l 2004) contains many OB and 
WN type stars, while the more evolved Quintuplet cluster and 
Central cluster with a ges of about 4 and 6 Ma (Figer et al.ll999al 
iPaumard et alll2006h contain WC stars as well. In addition, one 
well known LBV, the Pistol star ('Figer et al."l999a''^, and fur- 
ther LBV candidates (Geballe et al. 2000, Barniske et al. 2008) 
are found in the vicinity of the Quintuplet cluster. 

Massive stars in general, and WR stars in particular, are im- 
portant sources of ionizing photons, momentum and chemical 
elements to impact on the circumstellar medium. Stars that dis- 
play CNO-processed matter in a strong stellar wind are classified 
as WR stars of the nitrogen sequence (WN type). The cooler, 
late WN subtypes (WNL) usually contain some rest of hydro- 
gen in their atmosp heres, while the hotte r, early subtypes (WNE) 
are hydrogen free dHamann et al ] |l99lh . Typically, WNL stars 



are si gnificantly more luminous than WNE stars dHamann et al.l 
l2006l) . The WN stage may be followed by the WC stage, when 
the products of helium burning appear in the stellar atmosphere. 
Some evolutionary scenarios sugg est that the most m assive of 
these stars undergo LBV phases dLanger et al.l 199 4*). and the 
WNL phase might even precede the LBV stage dCrow ther 2003, 
ISmith & C onti 2008), whereas some evolutionary paths for stars 
with ini tial masses > 90 Mq seem to skip the LBV stage com- 
pletely (iMaeder et al.l l2008). However, the details of the evolu- 
tion of massive stars are still under debate. 

By means of quantitative spectral analysis in combination 
with comprehensive stellar model atmospheres the stellar param- 
eters can be assessed. This allows to determine feedback param- 
eters of the population of massive stars in the Galactic center and 
will provide constrains on stellar evolutionary models. 

However, the high extinction in this direction (Ay » 28 mag) 
prohibits observations in the UV and optical range. Therefore, 
IR spectroscopy is the prime tool to study massive stars in 
the GC. Different ESO IR-instruments at La Silla and VLT 



Paran a l were used to observe the Centra l cluste r dNaiarro et al 
1994', 'Krabbe et al ' 
2006, Martins et al 



19951 lEckart et al.1 |2004', 'Paum ard et al 
I2007h. the Arches cluster (Blum et a l.ll2001l 



Ma rtins et al and the Quintuplet cluster (iLiermann et al.l 

20091k 



* Based on observations collected at the ESO VLT (program 077.D- 
0281(A)). 



iNaiarro et all (1 1997b and iMartins et all (l2007l) both appl ied 
CMFGEN stellar atmosphere models dHilUer & Milled! 19981) to 
analyze /T-band spectra of massive stars in the Central cluster. 
They found that the stellar parameters of these stars are simi- 
lar to other Galactic WNL stars. Most of the WNL stars stud- 
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ied bv lMartins 6110(120071 note that some are classified only as 
Ofpe/WN9 stars) appear relatively rich in hydrogen. Therefore, 
the authors support the evolutionary sequence (Ofpe/WN9 <-> 
LBV) WN8 WN/C for most of the observed stars. 

In their s t udy of the brightest stars in the Arches cluster, 
iMartins et al.l (|2008|) found that these stars are either H-rich 
WN7-9 stars or O supergiants with an age of 2 to 4 Ma. The 
WN7-9h stars reveal high luminosities, log L/L© - 6.3, which 
is consistent with initial stellar masses of ~120Mo. The chem- 
ical composition shows both N enhancement and C depletion 
along with a still high amount of H, which leads the authors to 
conclude that the stars are core H-burning objects. It was found 
that the properties of the Ar ches massive stars argu e in favor of 
the evolutionary scenario of lCrowtheretaP (119951) O ^ Of ^ 
WNL + abs — » WN7. Thus, the massive stars in the Central and 
in the Arches cluster seem to conform to the standard evolution- 
ary models, and display stellar feedback parameters that do not 
significantly differ from other massive stars in the Galaxy. 

The present paper is based on our observations with the 
ESO-VLT Spectrograph for INtegral Field Observation in the 
Near-Infrared (SINFONI) of the central parts of the Quintuplet 
cluster. A spectral catalog of the point sources (hereafter LHO 
catalog) was presented bv Liermann et al.l(l2009h . The LHO cat- 
alog lists 13 WR stars, among them 4 WN and 9 WC types. In 
this paper we concentrate on the analysis of the WN stars, while 
the analysis of WC stars will be subject of a subsequent paper. 
Additionally, we analyze the star LHO 110, which was classified 
as Of/WN candidate in the LHO catalog. However, the spectrum 
of this star strongly resembles spectra of typical WN stars, there- 
fore we include LHO 1 10 in our sample. 

In order to analyze our program stars (WR 102d, WR 102i, 
WR102hb, WR102ea, and LHO 110), we fit their spec- 
tra with the Pots d am W olf-Rayet model atmospheres 
dHamann & GrafeneJ l2004l see also references therein) 
and derive the fundamental stellar parameters. By quantitative 
comparison of stellar ev olutionary models including rotation 
(iMevnet & Maedei1l2003h with the stellar parameters the initial 
masses and ages of the sample stars are determined. 

The paper is organized as follows. In Sect.|2]we give a short 
characterization of the applied model atmospheres. The spectral 
analyses is described in Sect.[3]and the derived stellar parameters 
are discussed in Sect.|4] Therein, mass-loss rates from radio free- 
free emission are debated (Sect. l4.2T l and our empirical results 
are compared with evolutionary models (Sect. l4.4l i. We conclude 
with a summary of the obtained results in Sect.|5] 

2. The models 

The Potsdam Wolf-Raye t non-LTE model atmospheres 
(iHamann & Grafene3l2004 PoWR) are based on the assump- 
tions of spherical symmetry and stationary mass loss. A velocity 
field is pre-scribed by the standard y6-law for the supersonic 
part, with the terminal velocity Uco as a free parameter The 
exponent /3 is set to unity for all models. In the subsonic region 
the velocity field is defined such that a hydrostatic density strat- 
ification is approached. Our models have the inner boundary 
set to a Rosseland optical depth of 20, thus defining a "stellar 
radius" /?» that represents the hydrostatic core of the stars. The 
"stellar temperature" r» is given by the luminosity L and the 
stellar radius via the Stefan-Boltzmann law, i.e. denotes 
the effective temperature referring to the radius R,. 

We account approximately for wind inhomogeneities 
("clumping"), assuming that optically thin clumps fill a volume 
fraction fy while the interclump space is void. Thus, the matter 



density in the clumps is higher by a factor D = /y ' compared to 
an un-clumped model of same parameters. All models are cal- 
culated with D-4 which is a rather conservative choice. Small 
scale random motions are considered by applying a microturbu- 
lence velocity of uu - 100 km s"'. Model atoms comprise hy- 
drogen, helium, CNO and iron-group elements. The latter are 
accounted for in a "superlevel " approach enabling blanketing by 
milhons of fines (iGrafener et al...2002i) . 

3. The analysis 

3.1. The observations 

We use flux-calibrated /iT-band spectra from the LHO catalog for 
the stars WR 102d (LHO 158), WR 102i (LHO 99), WR 102hb 
(LHO 67), WR102ea (LHO 71) and LHO 110. The data had 
been flux-calibrated by means of standards stars and then ap- 
plied for synthetic aperture photometry, yielding a mean error 
on the obtained K-band magnitudes of 0.4 mag (see LHO cat- 
alog for further details). For the comparison of the stellar pa- 
rameters, we refer t o previous studies of Galactic WN stars 
dHamann et al. 2006) and of the WN stars in the Arches clus- 
ter (Martins et al. 2008). Two further WR stars in the vicinity 
of the Quin tuplet cluster, WR 10 2c and WR 102ka, have been 
analyzed bv lBarniske et al.l (l2008l) recently. 

Our analysis is based on pronounced lines in the observed 
/iT-band spectra, which are due to helium and hydrogen: He i at 
2.059;um and 2.113;um, Hen at 2.189;um (10-7) , and Bry at 
2. 166 jjm blended with He ii (14-8) and various He i lines. A se- 
lection of /T-band lines is compiled in Table[T] 

3.2. PoWR - WN model fitting 

The normaUzed emission line spectra of Wolf-Rayet stars de- 
pend predominantly on two parameters, the stellar temperature 
Tt and a combination of parameters which was put in the form 
of the so-called "transformed radius" (ISchmutz et al.|[T989h 



' 'L2500kms-i / 10-4Moa-iJ ' 

The first step in the analysis is to adjust these two parameters 
such that the observed helium lines are reproduced. Concerning 
Bry, the hydrogen abundance is the relevant parameter, while the 
terminal velocity 0^ sets the width of the line profiles. 

Our analyses start from the two available grid^ of WN- 
type model atmospheres, calculated for a hydrogen abundance 
of 20% by mass ("WNL" grid) and for hydrogen-free stars 
("WNE" grid), respectively. Since all our program stars are of 
late spectral type with signatures of hydrogen in their spec- 
tra (WN9h), the WNL grid turns out to be more appropriate. 
This grid had been calculated for a fixed terminal wind veloc- 
ity of Doo - 1000km s After the best-fitting grid models were 
pointed out, we calculated individual models for each star with 
specially adjusted Uco and hydrogen abundance in order to opti- 
mize the fits. 

The results are shown in Table|2] while the observed spec- 
tra with the model fits are shown in Fig. [T] and Fig.|2] Note that 
we normalized the observed spectra by division through the red- 
dened model continuum. 

In the next step, we fit the spectral energy distribution (SED) 
which for our program stars is only available for the near- 
infrared range. Figure[T]shows the photometric fluxes in J, H and 

' http : //www ■ astro . physTk . uni -pot sdam . de/PoWR . html | 
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Table 2. Stellar parameters for the analyzed WN stars in the Quintuplet cluster 



LHO 
No. 


Alias 


Spectral 
type 


^s 
[mag] 


T, 
[kK] 


logR, 
[Re] 


[km s-'] 


Xh 
[%] 


[mag] 


Mk 
[mag] 


R, 
[Re] 


logM 

[Mo a-'] 


logL 

ILe] 


L/C 


67 


WR 102hb 


WN9h 


9.6 


25.1 


1.55 


400 


19 


8.0 


-8.3 


86 


-4.52 


6.42 


0.2 


71 


WR 102ea 


WN9h 


8.8 


25.1 


1.52 


300 


25 


6.1 


-8.3 


83 


-4.62 


6.39 


0.1 


99 


WR 102i 


WN9h 


10.1 


31.6 


1.65 


900 


45 


6.3 


-7.1 


42 


-4.79 


6.19 


0.5 


158 


WR 102d 


WN9h 


10.5 


35.1 


1.11 


700 


15 


7.1 


-7.0 


29 


-4.32 


6.07 


1.4 


110 




WN9h 


10.6 


25.1 


1.68 


300 


5 


8.7 


-7.6 


67 


-5.01 


6.20 


0.1 



Table 1. A selection of lines in the K-hand. 



^vac [jCm] 


Transition 


2.0379 


He II 1 5-8 


2.0431 


He I 6n 'P-4s ''S 


2.0587 


He I 2n ' P-2s ' S 


2.0607 


Hei 7d^D-4p^P 


9 0706/0809/084.9 


IV u — jij r 


2.1081 


Cm "in 'P 'S 
V- III jp r — j& o 


9 1 1 9^1 

jL. i IZU 


Up, 4.1; 3c Q„3p 
nc 1 H-s o — jp r 


2.1138 


Hei 4s 'S-3p 'P 


2.14999 


Hei7s3S^p^P 


2.1586 


Hei7p'P-4diD 


Bry blend: 




2.1614 


Hei7f3F-^d3D 


2.1623 


Hei7f 'F-^d'D 


2.1647 


Hei7g3G-4f3F 


2.1647 


Hei7g'G-4f 'F 


2.1652 


Hen 14-8 


2.1653 


Hei7d'D^f 'F 


2.1655 


Hei7d3D-4f3F 


2.1661 


Hi 7^ 


2.1792 


He II 23-9 


2.1821 


Hei7p3p-4d3D 


2.1846 


Hei7d'D~4p'P 


2.1891 


Hen 10-7 


2.2165 


He II 22-9 


2.2291 


Hei7s'S-^p'P 


2.2471/2513 


Niii5p2p-5s2s 


2.2608 


Hen 21-9 


2.30697 


Hei6p'P-4s'S 


2.3142 


He II 20-9 


2.3470 


Hen 13-8 



is: from 2MASS dSkrutskie et al.ll2006h . and our flux-calibrated 
spectra from the LHO catalog together with the photometric 
magnitude derived there. Crowding in the cluster field does not 
allow photometric fluxes to be extracted, e.g. from Spitzer IRAC 
or MSX images. 

In the case of WR 102d we encountered a problem with 
the 2MASS magnitudes from the final release of the catalog, 
which gives unplausible values. Moreover, we find a coordi- 
nate offset between the 2MASS position and the LHO posi- 
tion of WR 102d of almost 5 "; we identify the 2MASS posi- 
tion rather with LHO 147 (see LHO catalog). However, the data 
from the 2MASS intermediate release catalog, although consid- 
ered to be obsolete, contain a point source at exactly the position 
of WR 102d with fluxes that fit very well to our calibrated K- 
band spectrum. Therefore we adopt these values for the SED 
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Fig. 1. Spectral energy distributions for the sample stars. Blocks 
with labels indicate the 2MASS magnitudes for the J-, H- and 
/Ts-band (dark blue boxes), as well as the LHO magnitude 
(light blue box). The solid line is the flux-calibrated LHO spec- 
trum while the dotted line refers to the reddened model contin- 
uum. 



fit. The comparison with the near-IR photometry presented by 
Figer et al. ( 1999a, filters J, H2 and K2) supports this choice. 

The model continuum (red dotted lines in Fig.[TJ is then fitted 
to the observations by adjusting two parameters, the reddening 
parameter (e. g. given i n the form of the color excess Eb-v in the 
narrow-band system of Smith (1968), see Table|2]i and the stellar 
luminosity. We adopt a distance to the Quintuplet of 8 kpc, cor- 
responding to a distance modulus of 1 4.52 mag ( Reid 1993). The 
reddening law which we adopt is from lMoneti et al.i(i200L) . The 
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reddening parameter translates to the usual Johnson sys- 
tem as Eb-v - 1-21 X Eh-v, while the extinction in the /T-band 
is Ak - 0.42 X Eh-v Thanks to the approximate invariance of 
Wolf-Rayet line spectra for models with same transformed ra- 
dius, models can be scaled to different luminosities for the SED 
fit while the line fit is preserved. Generally we can reproduce 
the SED quite well. We attribute the remaining discrepancies to 
problems of the 2MASS photometry in the crowded field. 

Remember that we have already employed the model con- 
tinuum to normalize the observed spectra. Insofar, the described 
two steps of the analysis are in fact an iterative procedure. The 
derived extinction parameters are compiled in Table|2] and fur- 
ther discussed in Sect. 13. 71 After we fitted the line spectrum and 
the SED, the whole set of stellar parameters is established (see 
Table|2ll. 

3.3. WR 102ea and WR 102hb 

These stars, as well as LHOllO (see Sect. l3.5l l. do not show 
clear He n lines in their spectra, which is indicative for a low- 
temperature regime in their atmospheres. Attempts to fit the 
spectra with models of too cool temperatures makes those lines 
disappear completely and overestimates the He i lines instead, 
especially Hei 2.059yum. In comparison, a hotter model pro- 
duces too strong He ii lines in a strength in which they are not 
observed. 

Additionally, these stars have terminal wind velocities of 
300-400kms which is rather slow for WR winds. In com- 
parison, the adopted standard microturbulence - 100 km s ' 
might be too high. Therefore, we tested different microturbulent 
velocities down to uu = 20kms"'. However, the best fit to the 
observed profile shapes is obtained for uu - 80kms In any 
case, the effect of this parameter on the synthetic /T-band spectra 
is small. 

3.4. WR 102dand WR 102i 

The spectra of these two stars both show Hei and Hen lines 
which allow the temperature to be determined more reliably 
than in the previous cases. As for WR 102i, however, our model 
fall short in reproducing the strong emission feature of He i at 
2.113jum, while the absorption feature is too strong. Although 
a bit hidden in the telluric residuals, the same seems to hold 
for He: at 2.059/im, but the Hen absorption at 2.189yum fits 
perfectly. Alternatively, we found slightly hotter models which 
reproduce the Hei emission, but then fail with the Hen lines 
turning into emission. We are not sure about the reason for 
these inconsistencies, but we note that e. g. the N iii doublet at 
2.247/5 lyum seems to show a different radial velocity, possibly 
indicating that the spectrum is composite from a binary. This 
needs further observations to be confirmed. For the moment, the 
fit problems introduce extra uncertainty which might be lower 
for the alternative models. 

In the case of WR 102d the lines He i at 2.059 fim and He ii 
at 2.189//m can be reproduced consistently by the model. With 
the stellar temperature being well constrained, the Bry fit also 
yields an accurate hydrogen abundance (see Fig.O. 

However, the spectra of both stars show remarkably strong 
Niii lines at 2.247/5 ljum. Although we already spent some ef- 
forts in improving the corresponding part in our nitrogen model 
atom, our models notoriously fail to reproduces this doublet in 
the strength it is observed here, while it is reasonably well fitted 
for the weak emission in WR 102ea and WR 102hb. 



In comparison to the other stars in the sample these two stars 
show the highest terminal wind velocities, which are more rep- 
resentative for WN stars. The final best fitting models were cal- 
culated with the standard microturbulence i>d - lOOkms 

3.5. LHO110 

LHO 110 was classified as 06-8 Ife in the LHO catalog. The 
emission line spectrum is very similar to the spectral appear- 
ance of the WN stars, although the spectral features are less 
pronounced. In detail, the star has a similar Ks magnitude as 
the WN9 stars and also shows similar line ratios for (Hen 
2.189A/m)/(Hen/Brr 2.166 jum) and (Br7)/(Hei 2.059//m). 

We analyzed the star in the same way as described above 
for the WN stars. We calculated tailored models to fit the line 
spectrum and SED accordingly, see Fig. [T] and |2] The obvious 
weakness of Bry leads to a hydrogen mass fraction as low as 
5%. However, models with zero hydrogen cannot fit the blend at 
all which excludes the star to be hydrogen free. Hence the spec- 
tral classification in the LHO catalog must be revised to W N9h. 
Following the naming conventions of van der Huchtl (1200 1 ) this 
star could become WR 102df. The flux-caUbrated spectrum can 
be well fitted with the model SED, but cannot be connected to 
the 2MASS photometry. As in the case of WR 102d, we attribute 
this to the problems of the 2MASS measurements. The star lies 
in the central region of the cluster where crowding might have 
affected the 2MASS data. Further on, a set of small absorption 
lines is present in the spectrum, which we cannot identify. They 
might indicate that LHO 1 10 is a binary. However, like the other 
stars in our samp le, LHO 110 has not been detected in X-rays 
(Wang et al."2006'). Collision of stellar winds in binary WN stars 
often make them brighter X-ray sources than single stars (e. g. 
lOskinova 2005). Therefore, the binary nature of LHO 110 can- 
not yet be confirmed. 

3.6. Error margins 

The temperature of the model to be fitted can quite accurately be 
determined for at least those sample stars that show Hei and 
(even weak) Hen lines. In comparison, too cool or too hot a 
model would produce He n lines in a weakness or strength in 
which they are not observed. From this, we consider our temper- 
atures to be exact within the range of ±0.025 dex, in general. 

Different influences contribute to the error of the luminosi- 
ties given in Table|2] The fit of the SED depends on the accuracy 
of the photometric measurements. The 2MASS colors often suf- 
fer from source confusion, while the magnitudes from the 
LHO catalog have been calibrated to an accuracy of +0.1 mag. 
The reddening parameter Eb-u obtained from the slope of the 
SED is accurate by about +0.1 mag, which gives only a 0.04mag 
uncertainty in the /T-band extinction. 

Given that we fit the SED in the Rayleigh- Jeans domain, the 
bolometric correction factor (10** *^'') scales with tI. Hence the 
temperature contributes about 10% to the error margin of the 
luminosity. In combination we think that log L has a typical error 
of about ±0.16 dex. 

As described in Sect. 13. 21 the hydrogen content was deter- 
mined by fitting models with different mass fractions of hydro- 
gen to the spectra. Fig. |3] shows the Bry/ Hei/ Hen blend at 
2.166/im of WR 102d. The model with a hydrogen abundance 
of 15 % by mass gives the best fit to the observation (dotted line) 
within the error range of 10 to 20% hydrogen. Therefore we 
conclude that the derived values for are correct within ±5 % 
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Fig. 3. Normalized line spectrum of WR 102d showing the 
Bry / He i/ He ii blend (blue solid line) with the best fitting model 
(red dotted line). For comparison we plot models with 20% 
(green dash-dotted line) and 10% (black dashed line) hydrogen 
by mass but otherwise identical parameters. 



by mass, accounting also for uncertainties in the fit of the con- 
tinuum. 



3.7. Extinction 

In the LH O catalog \ye assu med an average extinction of A^: = 
3.28 mag (iFiger et alJl999ai) for the Quintuplet cluster. With the 
model fit of the SEDs for our sample stars, we have determined 
individual values for the reddening. The results for the stars are 
listed in Table|3] with an aver age of A;^ = 3.1 + 0.5 mag, sim- 
ilar to the one determined bv iBarniske et al. (2008;), who ana - 
lyzed the stars WR102c and WR102ka, and .Figer et aPdlQQQal) . 
The average extin ction of the Arches cluster was found to be 
Ak - 2.8 mag bv iMartins et all (l2008l) . shghtly lower than for 
the Quintuplet. 

Within the errors, we conclude that the major part of the ex- 
tinction can be attributed to the foreground interstellar redden- 
ing. From the limited number of stars analyzed it is not possi- 
ble to derive an extinction map for the Quintuplet cluster. We 
find a slight increase in the determined values towards LHO 1 10 
which is located closer in the cluster center. This could be due 
to the intrinsic cluster extinction expected from the presence of 
dust-producing WC stars (.Tuthill et al...2006.) . 



4. Discussion 

4.1. Steiiar parameters 

Fig.|4] shows the derived empirical mass-loss rates versus stel- 
lar luminosity. The stars analyzed in this work fall below the 
mass loss-luminosity relation of lNugis~& Lamers (2000) for 
stars containing 40 % hydrogen per mass (lower line). From their 
determined hydrogen abundance (see Table|2]l, one would expect 



Table 3. Extinction derive d from fitting the m odel SEDs, ap- 
plying the reddening law of iMoneti et alJ (1200 ih in the infrared 
range. 



Object 


Av 


Ak 




ref. 




[mag] 


[mag 






WR102hb, Q8 


30 


3.4 




this work (t.w.) 


WR102ea, Q 10 


23 


2.6 




t.w. 


WR102i 


24 


2.7 




t.w. 


WR102d 


27 


3.0 




t.w. 


LHO 110 


33 


3.7 




t.w. 


mean 


27 ±4 


3.1 ± 


0.5 


t.w. 


WR102C 


26 ± 1 


2.9 ± 


0.1 


Bamiske et al. (20081 


WR102ka 


27 ±5 


3.0 ± 


0.6 


Barniske et al. (2008) 


Quintuplet mean 


29 ±5 


3.28 


t 0.5 


Fiser et al. (1999a) 


Arches mean 


24.9 


2.8 




Martins et al. (2008) 
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Fig. 4. Mass-loss rate versus stellar luminosity for the sam- 
ple stars. For comp arison we show Arches WNL stars from 
IMartins et all (l2008l labeled with F a nd B numbers), WR102c 
and WR102ka fro m Barnis ke et al] (12008 ) and Galactic WN 
stars dHamann et al.i 2006. open symbo ls). The solid lines refer 
to the mass loss-luminosity relation of iNugis & LamersI (l2000h 
for hydrogen-free stars (green upper line) and stars containing 
40 % hydrogen per mass (red lower line). 



to find these stars positioned between this relation and the one 
for hydrogen-free stars (upper line). 

Not only our program stars have rather low mass-loss 
rates compared to their luminosities, the same holds for 
other Galactic WN stars shown in Fig.|4] (open symbols, 
from iHamann et al. I2006h. i ncludi ng the WNL stars in the 
Arches cluster ( Martins et al.l l20()8l symbols labeled B and F, 
from) which have similar hydrogen contents. For radiation- 
driven winds , the mass-loss rates scale with the metallicity 
dGrafener & H amann 2008). Therefore, the relatively low mass- 
loss rates provide an indirect argument that the metallicity in our 
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Fig. 5. HRD showing the analyzed WN9h stars (red filled sym- 
bols). For compariso n we plot single Galactic WN stars from 
iHamann et alJj 20061) (open symbols), WNL stars of the Arches 
cluster from iMartins et al.l (2008, labeled with F and B num - 
bers), and WR102c and WR102ka from Barni ske et al] (|2008|) . 
Different symbols refer to different spectral subtypes, cf. inset. 
Hydrogen and helium zero-age main-sequence are given as solid 
lines. 



program stars is not significantly higher than in other Galactic 
WN stars studied so far 

However, the possibility that the Galactic center migh t be 
metal enriched is still under debate. iNaiarro et al.l (|2004|) ob- 
tained solar metallicities for a sample of stars in the Arches clus- 
ter, while recent results on a larger sample in the same cluster 
favored slight metal enrichment in the range of 1.3 to 1.4Zq 
(Martins et al."2008'). For the Quintuplet cluster, 'Najarro et alj 
(2009) used quantitative spectral analysis of the LBV Pistol 
star and LBV candidate FMM362 to measure the metallicity 
in the cluster Their results indicate solar iron abundance and 
roughly twice the solar abundance in the a-elements. They sug- 
gest that the enrichment in a-elements versus Fe found in the 
two Quintuplet LBVs is c onsistent with a top-h eavy initial mass 
function (IMF) in the GC. iDavies et all (l2009h analyzed the red 
supergiant GMM7 (LH07) in the cluster The authors find a 
slight iron enrichment which they explain by the evolved nature 
of the object and the resulting hydrogen depletion. They argue 
that the "initial" iron abundance would be in agreement with so- 
lar values. 

The Hertzsprung-Russell diagram (HRD, Fig.|5]i shows the 
program stars of the Quintuplet cluster in comparison to other 
Galactic WN stars. W e include WN stars from the Arches cluster 
(iMartins et alJl2008l symbols labels B and F, from) and the two 
bright stars WR102c and WR102ka (iBarni ske et al. 200j). The 



whole sample of single Galactic WN stars (iHamann et aLll2()()a, 
open symbols) is also shown for comparison. 

Obviously, the GC stars populate a special domain in the 
HRD. They lie within a group of rather cool, T, < 50 000 K, 



but very luminous WN stars with log (L/Lq) > 6.0. Like most 
WNL stars in the Galaxy, our GC stars show some rest of hydro- 
gen with mass fractions of the order of 20%. While about half 
of the Galactic WN stars belong to the "early" subclass (WNE) 
that resides to the left side of the ZAMS (see Fig.|5]l, no single 
such hot WN star has been found in the Galactic center region. 

One may wonder if the low stellar temperatures obtained for 
the GC stars are due to systematic errors of spectral analyses that 
are based o n the A"-band only. Bu t the stars in the Arches cluster 
analyzed bv I Martins et alJ (l2008h show that independent groups 
arrive at similar results. 

Can the lack of WNE stars be attributed to selection effects? 
Compared to WNL stars, WNE stars are roughly th ree magni- 
tudes fainter in the A'-band (cf. IHamann et alj|20()6l the differ- 
ence is similar as for the visual band). Since the apparent K- 
band magnitude of our program stars is about 10 mag (Table|2]i, 
the WNE stars should not all have escaped from detection in the 
LHO catalog which is complete to the 13* magnitude. 

A statistical comparison with the WNE star population in 
other regions of the Milky Way confirms that WNEs are usu- 
ally detected with slightly fainter /T-band magnitudes compared 
to their local WNL counterparts. However, the number ratio 
■ ^WNL : -^WNE var ies largely from 24:33 for Galactic WN 
(IHamann et alJ2006l) . 4:1 1 fo rWd 1 jCr owther et al. 2006), 10:3 
for the GC Central cluster dMartinset all 12007 ). 13:0 for the 
Arches (Martins et al. 2008), and 5:0 for the Quintuplet, indi- 
cating a lack of WNE stars in the GC region. 

Hence, we must conclude that the population of WN stars 
in the Galactic center region is different from the Galactic WN 
population in general. This may be attributed to various reasons: 
a special star-formation history, a different initial mass func- 
tion (IMF), or a (slightly) higher metallicity that increases the 
mass loss during stellar evolution. A closer investigation of this 
question may give important insight into the conditions of the 
Galactic center environment. 

For all stars in the sample we derive the wind momentum 
efficiency number 77 - Mnoac/L, see Table|2l Recall that we have 
adopted a clumping contrast of D = 4 throughout this paper (cf. 
Sect.|2]i. If D is higher, this would reduce the empirical mass- 
loss rates oc D ''^ and thus lead to smaller 77 values. However, 
this number cannot exceed unity in r adiation-driven winds unless 
mult iple scattering is taken account. [Grafener & Hamarnil (I2005L 
I2OO8 ) have shown that 77-values of a few can be easily obtained 
in hydrodynamically consistent WR models with full radiative 
transfer. For all stars we find values of 77 close to unity, which is 
consistent with these stellar winds being radiatively driven. 

4.2. Radio mass-loss rates 

Three of the sam ple stars are known to be radio sources from 
iLang et al] (l2005h . We calculated the mass-loss rates from ra- 
dio fre e-free emission with rad io fluxes from these authors and 
foUow lWright & Barlowl d 19751) 



M = Do, 



(23.2)3/4 (ygffV) 



1/2 



(I)- 



(2) 



The free-free Gaun t facto r gff is derived from the relation by 
iLeitherer & Robert! (Il99ll) with an assumed electron temper- 
ature of lOOOOK. Terminal velocity Uoo and mean molecular 
weight jj are taken from our analysis (see Table|4]i. Assuming 
that helium stays ionized in the radio emitting region of the stel- 
lar wind, the mean number of electrons is set to y = 1 . Radio 
mass-loss rates were calculated from the radio flux at 22.5 GHz 
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Table 4. Mass-loss rates derived from radio free-free emission 
(Miadio), compared to those determined from the ^T-band spec- 
tra i n the presen t paper (Mir). Based on radio measurements 
from iLang et aP 120051) . we derive Mjadio consistently with the 
parameters determined in this paper. The values in parenthe- 
sis are the mass-loss rates from lLang et alj (l2005h who adopted 
Ilea - 1000km s"' and fJ. -2 for all stars. 



Object 


alias 




A* 


log M,„d,o 


log Mir 






[km s-'] 




[Mo a-'] 


[Mo a-'] 


WR102ea 


QR5 


300. 


2.30 


-4.44 (-3.96) 


-4.62 


WR102d 


QR8 


700. 


2.78 


-4.30 (-4.37) 


-4.32 


LHOllO 


QR4 


300. 


3.52 


-4.12 (-3.82) 


-5.01 



for QR5 (WR102ea) and QR8 (WR102d), and at 8.5 GHz for 
QR4 (LHO 110). This allows the derive d mass-loss rates t o be 
compared directly with those obtained bv lLang et al.l (l2005h . but 
note that their original values for M, given in parenthesis in 
TableH] were calculated with different assumptions on the ter- 
minal velocity and the mean molecular weight. 

For QR5 (WR102ea) and QR8 (WR102d) the mass-loss 
rates which we deduce from the radio measurements by 
iLang et al. (2005') are in excellent agreement with those obtained 
in the present paper from the ^T-band analysis (see Table|4|. 
Thus, these two radio detections can be confirmed as stellar wind 
sources. 

It has to be noted that the radio mass-loss rates quoted in 
Table|4] are derived for an unclumped medium (clumping con- 
trast D - I), while the model TiT-band spectra are calculated for 
slightly clumped winds (D - 4). Hence, the good agreement 
implies that the degree of clumping generally decreases from 
the line forming region towards the radio emitting a s previously 
found for OB stars (Puis et al. 200^ and WR stars (iNugis et al.l 
ll998llLiermann & Hamann.2008iK 

In the case of LHO 110, however, the mea sured radio flux is 
about ten times higher then expected from the Wright & Barlow 
(1975) model compared to the mass-loss rate derived from the 
/T-band analysis. This might indicate that the radio emission 
is produced in a colliding-wind zone. On the other hand, bi- 
naries often show a non-thermal radio spectrum with a nega- 
tive radio spectral index. For LHO 1 lO. ILang et al.l (|2005) found 
a s pectral index of -1-1.4, w hich is closer to the prediction of 
the IWright & Barlowl ( Il975h theory for thermal wind emission 
(spectral index -1-0.6). Considering these contradictory indica- 
tions, the nature of the radio emission and the binary status of 
LHO 110 remains unsettled. 

4.3. LBV candidates? 

We find two stars in our sample, WR102ea and WR102hb with 
very high luminosities and low terminal velocities. Their mass- 
loss rates are comparable to those derived for the Pistol star and 
FMM362 (Najarroetal. 2009), or WR102ka dBarniske et alj 
I2OO8 ) thus leading to the conclusion the stars might be LBV 
candidate stars. For LHOllO the slightly smaller luminosity 
sets the star in the HRD just below the Humphreys-Davidson 
limit adapted from'Figer et al.l (Il998i their Fig 12), while other 
Quintuplet LBV candidates are above this empirical relation (see 
Fig.|6]l. Also the very low mass-loss rate might question the LBV 
scenario for LHO 1 10, despite its spectral appearance. 



; ] ] ^ ] ] ] ^ ] ] r 




iog(r./K) 

Fig. 6. HRD with LBV (candidate) stars in the vicinity of the 
Quintuplet cluster In addition to the program stars of this paper 
we indicate the Pist ol star and FMM 362 (parameters taken from 
iNaiarro et al.ll2009h and WR 102ka dBarniske et al J l2008l). The 
Humphreys-Davidson limit (solid line, adapted from'Fig er et alJ 
1998 D is shown a s wel l as the hot LBV minimum light strip 
from IClark et al.1 (l2005l dashed, extrapolated with thin dashes 
to higher l uminosities ) . Eve n the steeper version of this latter 
Umit from lGroh eTaP d2009L dash-dotted line) puts WR 102ea, 
WR 102hb, LHO 1 10, and WR102ka to its right side, i.e. in the 
suggested LBV domain. 

'Clar k et aP d2005h proposed an empirical limit for LBV 
(candidate) stars in the quiescent state. This hot LBV minimum 
light strip is defined up to luminosities of log(L/Lo) ~ 6.2 by 
observed LBVs and LBV candidates, but a possibl e extension 
towards higher luminosities is discussed. Further on. lClark et al.l 
(l2005h argue that LBV stars in this part of the HRD will look 
like WNL stars. As recently presented by'Groh et al.' ('2009'), this 
minimum light strip is connected to the strong-variable LBVs 
suspected to be fast rotators reaching their critical rotational ve- 
locity, e. g. the known strong variable LBVs AG Car and HR Car 
lie at this limit. In contrast slow-rotating LBV stars are assumed 
to be less variable. In any case, we don't find hints in the spectra 
of our sample stars to be fast rotating and the question of vari- 
ability cannot be settled without further observations. Still, all 
three stars, WR 102ea, WR 102hb, and LHO 110, are thus left in 
the possible LBV regime within the HRD. 

WR102ea and WR102hb both are placed above the evolu- 
tionary track with 120 initial mass, see Sect.l4~4l and Fig.|7] 
But as described recently by Ma eder et al.l (l2008h . such massive 
stars might skip the LBV phase in their evolution completely. 

4.4. Stellar evolution 

The H RD with evolutionary tracks from iMevnet & MaedeJ 
d2003h is shown in Fig.|2l Most sample stars lie within the tracks 
for solar metallicity between 60 and I2OM0 initial mass imply- 
ing that the WN stars are descendants from the most massive O 
stars. The tracks without rotation (dashed lines) extend further 
into the cooler temperature domain than the tracks with rotation 
(solid lines). This might imply that the initial equatorial rota- 
tional velocity of the sample stars was smaller than the assumed 
300km s"' in the evolutionary models. Similar results were al- 
ready encountered when comparing the Galactic WN stars with 
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Fig. 7. HRD with evolutionary tracks from iMevnet & Maederl 
(l2003h for Mini of 60, 85 and 120 Mq. The dashed hnes re- 
fer to tracks without rotation, soHd Hnes to tracks with t^ini - 
300km s~'. Symbols are for the analyzed Quintuplet WNL stars 
as in Fig.|5] 

the Geneva evolutionary tracks (iHamann et al ] l2006h . This indi- 
cates that the stellar evolution of massive stars is not yet under- 
stood fully. Two stars in the sample, WR 102hb and WR 102ea, 
lie above the most massive track with 120 Mq initial mass. We 
compare the stellar pa r amete rs of the stars with evolii t ionary 
tracks by iLanger etal] d 19941) as shown by .Figer et alj (I1998L 
Fig. 15a). Thus the initial stellar mass for these stars would be in 
the range 150 < Minit < 200 Mq. 

For those stars that can be assigned to Geneva evolutionary 
tracks, we estimated the age of the stars to be 2.1 to 3.6 Ma and 
find present-day masses of 30 to 45 Mq. This is in agreemen t 
with masses we derive following the relation bv lLangeij(ll989h . 
which is strictly for pure helium stars only. Considering the de- 
termined stellar temperatures and luminosities we find from the 
Geneva models that the stars are still hydrogen burning objects, 
in spite of their small hydrogen abundances. 

5. Summary 

We present the analysis of five WN9h stars in the Quintuplet 
cluster. ^T-spectra were taken from the LHO catalog and fitted 
with tailored Potsdam Wolf-Rayet models for expanding atmo- 
spheres (PoWR) to derive the fundamental stellar parameters. 
All studied stars are very luminous (log (L/Lq) > 6.0), while 
their stellar temperature corresponds to their late WN subtype 
(T, < 35 kK). Along with other WNL stars in the Milky Way 
the Quintuplet stars form a distinct group that occupies a spe- 
cific region in the HRD, separate from the WNE stars. They still 
contain a significant amount of hydrogen (up to 45 % per mass) 
in their atmospheres and show typical WN-like mass-loss rates. 

Radio detections exist for three of our sample stars. In two 
cases, WR 102d (QR8) and WR 102ea (QR5), the radio emis- 



sion agrees with the prediction for free-free emission from their 
stellar winds. In contrast, LHO 110 (QR4) is much brighter in 
the radio than consistent with these assumptions, which possibly 
indicates that it might be a binary with colliding winds, although 
none of the sample stars has been detected in X-rays. 

The two most luminous stars of our sample correspond to ini- 
tial masses above 150 M©. The other three, less luminous stars 
can be compared with detailed evolutionary tracks, revealing ini- 
tial masses in the range of 60 to 120 Mq, and ages of about 2.1- 
3.6 million years. As in the case of other Galactic WN stars, 
the evolutionary tracks without rotation seem to fit better to the 
stars' location in the empirical HRD. According to the evolution- 
ary models, these stars are still core-hydrogen burning objects. 

From the analysis of the individual stars we obtained the in- 
terstellar reddening and extinction. Based on the small sample, 
no extinction map could be established but an average extinc- 
tion of Ak - 3.1 ± 0.5 mag (Ay - 27 ± 4 mag) was determined 
towards the Quintuplet cluster. 
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